Non-local state-swapping of polar molecules in bilayers 
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The observation of significant dipolar effects in gases of ultra-cold polar molecules typically de- 
mands a strong external electric field to polarize the molecules. We show that even in the absence 
of a significant polarization, dipolar effects may play a crucial role in the physics of polar molecules 
in bilayers, provided that the molecules in each layer are initially prepared in a different rotational 
state. Then, inter-layer dipolar interactions result in a non-local swap of the rotational state be- 
tween molecules in different layers, even for weak applied electric fields. The inter-layer scattering 
due to the dipole-dipole interaction leads to a non-trivial dependence of the swapping rate on den- 
sity, temperature, inter-layer spacing, and population imbalance. For reactive molecules like KRb, 
chemical recombination immediately follows a non-local swap and dominates the losses even for tem- 
peratures well above quantum degeneracy, and could be hence observed under current experimental 
conditions. 

PACS numbers: 67.85.-d, 34.50.Cx 



A new generation of experiments has started to explore 
the remarkable novel physics of dipolar gases, in which 
dipole-dipole interactions play a key role [1| . These inter- 
actions, being long-range and anisotropic, differ from the 
short-range isotropic interactions which have dominated 
up to now ultra-cold atomic physics. Polar molecules 
are expected to provide fascinating new scenarios for 
quantum gases due to their large electric dipole mo- 
ments. Recent experiments on preparation and con- 
trol of ro-vibrational and hyperfine states of KRb at 
JILA [2j open new perspectives towards a degenerate 
quantum gas of polar molecules. Unfortunately chemical 
recombination due to the reactive character of KRb has 
up to now prevented to reach quantum degeneracy Q. 
However, dipolar interactions between partially polarized 
KRb molecules have been recently shown to significantly 
reduce chemical recombination in constrained quasi-2D 
geometries Q. 

Although a polar molecule in its lowest ro-vibrational 
state may have a large dipole moment in the molecu- 
lar frame, in the absence of an external electric field the 
dipole moment in the laboratory frame averages to zero. 
As a consequence, the observation of dipolar effects in 
gases of polar molecules typically demands the polar- 
ization of the molecules in an external electric field 
For example, a KRb molecule in its singlet ground state 
has a permanent dipole moment of d = 0.566 Debye in 
the molecular frame [2J , but rather large fields of several 
kV/cm must be employed to reach effective dipoles of 
^0.2 Debye in the lab frame [f|. 

However, remarkably, dipolar interactions may play a 
significant role even in the absence of an external electric 
field, if the molecules are prepared in different rotational 
states [6j . In that case, the interaction between molecules 
in different rotational states is produced by exchanging a 
quantum of angular momentum (state swap) , resembling 
the resonant interaction of an electronically excited atom 
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FIG. 1: Bilayer set-up with molecules in rotational state 
ipo.o (layer A) and ipi^ (layer B). A state swap has been 
sketched. 

and a ground state one. 

In this Letter, we show that bilayer systems may pro- 
vide an excellent environment for observing the non- 
trivial physics associated with dipole- induced state-swaps 
in current experiments. We consider polar molecules in 
a bilayer geometry, without inter-layer tunneling, pre- 
pared such that the molecules in layer A (B) are in 
the ground (first-excited) rotational state (Fig. [1]). The 
dipole-dipole interaction induces a molecule in one layer 
to interchange its rotational state with a molecule in a 
different state in the other layer, as sketched in Fig. [T] 
Such a process resembles spin-exchange in the collision of 
polarized atomic beams [7( and spin-changing collisions 
in spinor gases We show that the special features of 
the inter-layer scattering lead, especially for quantum de- 
generate samples, to a non-trivial dependence of the swap 
rate on temperature, molecular density, inter-layer spac- 
ing, and population imbalance. Interestingly, for reac- 
tive molecules like KRb chemical recombinations follow- 
ing a state-swap dominate losses even well inside the non- 
degenerate regime and could be observable under current 
experimental conditions. 

We are interested in fcrmionic polar molecules con- 
fined in two thin parallel layers, A and B, separated by 
a distance A along z, and with negligible hopping be- 
tween them (Fig. [1]). We consider diatomic molecules 
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with a permanent electric dipole moment d in the elec- 
tronic singlet state ( 1 S) and in the vibrational ground 
state. A weak external electric field E = £ z is applied. 
The rotational states of the molecules are described by 
the Hamiltonian 



H mo \ = BJ — d ■ E, 



(1) 



where B is the rotational constant, J is the angular mo- 
mentum operator, and d = d(sin 9 cos ip, sin 9 sin ip, cos 9) 
is the dipole moment operator, where the angles 9 and tp 
describe the orientation of the molecule in the laboratory 
frame. For a weak electric field, such that p = d£ /B -C 1, 
the eigenstates ipj.M of the Hamiltonian (JTJ are still well 
described by the quantum numbers J and M associated 
with the angular momentum J, and they can be found 
using perturbation theory [t| [l(| . The electric field lifts 
the degeneracy of the three J = 1 states, which become 
separated by an energy AE = Ei, — £i,±i = 3f3 2 B/20. 
The effective dipole moment in the laboratory frame, (d) , 
is very small for a weak field. For example, the effective 
dipole in the rotational ground state is (</?o,o|d|yo,o) — 
fid/Zi, giving 0.019 Debye for KRb and /?'= 0.1,' cor- 
responding to £ = 0.4kV/cm; this results in a dipolar 
interaction five times smaller than that for atomic mag- 
netic dipoles in Chromium jl|. 

We suppose that, initially, the molecules in layer A are 
prepared in the state <po,0: while those in B are prepared 
in ipifi, as illustrated in Fig. [TJ Such selective prepa- 
ration may be achieved by applying a spatially varying 
electric field, and then transferring the molecules of layer 
B into the desired state with a microwave transition [11| . 
The dipole-dipole interaction between two molecules with 
dipole di i2 at positions rx 2 is 



Hdd — 
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di-d 2 3(di-r)(d 2 -r) 



with r = ri— r 2 . Since the anisotropic dipole-dipole inter- 
action does not preserve the rotational quantum number, 
a collision may lead to an interchange of the rotational 
state of molecules in different layers: 



(3) 



The matrix elements of Hdd involving the states (^o,o 
and ipifl are given by, neglecting terms 0(f3 2 ), 



3r 3 



3z 2 



(4) 

V dd {v). (5) 



Transitions into rotational states other than <po,o and ip\,o 
can be neglected. They involve an energy change of AE, 
which is much larger than other energy scales in the sys- 
tem (the rotational constant B is typically very large), 
such as the dipolar energy Ed = e? 2 /(47reoA 3 ). For KRb, 



for example, B/h = 1.11 GHz @, thus AE/k B = 80^K 
is much larger than £x>/fce = 15nK for A = 532nm. 
Therefore, we only consider processes between the states 
Vo,o and <p lfi [14 1. 

At this point the dipolar interaction is easily diag- 
onalized using the symmetric and antisymmetric two- 
body states \S,A) = (\(p£ fi <p?, a ) ± Ko^o))/^. for 
which (A\H DD \A) = -V dd {r), (S\H DD \S) = V dd {r), and 
(S\Hdd\A) = 0. The inter-layer scattering is hence de- 
composed into two separate channels, namely the sym- 
metric one, with interaction potential V dd (r), and the an- 
tisymmetric one, with — V dd (r). It may be noted that in 
analogy with spin-exchange collisions between polarized 
atomic beams [7j , the interaction between two molecules 
can be written as V — V dd (r)Ps — V dd (r)PA, where 
Ps,Pa are the projection operators for states \S), \ A). 

We introduce the scattering amplitudes /+ and /_, 
which characterize the low-energy inter-layer scattering 
for the symmetric and antisymmetric channel, respec- 
tively. The scattering amplitude for the state-changing 
process in Eq. ^ is then given by 



/«(*) = £(/+(*:) -/-(*)), 



(6) 



with hk the relative momentum between the colliding 
particles. The amplitude for state-preserving collisions, 

<P0,0 VIA -> ^0,0 Vlfii is /bp = W- + /+)• 

Collisions involving particles in different layers can 
be described here by two-dimensional scattering (s-wave 
only). The amplitudes /+ and /_ are found by numer- 
ically solving the 2D Schrodinger equation with the po- 
tential ±V dd (x, y, z = A). This may be compared, using 
Ref. to the result from the Born approximation for 

f±, 



(2) |/ sc (g)| 2 = C/ V[2 + 7 r(L 1 (2g)-/ 1 (2 g ))] 2 + O(C/ 4 ), (7) 



with q = kX and Li(x) (Ii(x)) the modified Struvc 
(Bessel) function. We use the unit of energy Eq = 
H 2 /(m\ 2 ), with m being the mass of the molecule, and 
the dimensionless coupling strength Uq = Ed/(3E ). 
To be specific, we will take the values for KRb and 
A = 532nm for the plots, so Eo/ks = 13.5nK and 
Uq = 0.38. In Fig. HJa) we show that the expression in 
Eq. ([7]) is in good agreement with the numerical solution 
of the Schrodinger equation, for the regime of interest. 

Note that |/ sc | 2 presents a non-monotonical depen- 
dence with the relative kinetic energy, and falls off slowly 
for large energies, see Fig. (5Ja). Interestingly, the swap- 
ping amplitude in Eq. (JTJ) is 1/3 of the scattering ampli- 
tude expected for fully polarized molecules calculated in 
the same approximation. On the other hand, the scat- 
tering amplitude for state-preserving collisions is much 
smaller. Since there is no direct process (the matrix el- 
ement vanishes, Eq. (0|), these collisions are of higher 
order: \f sp \ 2 = (D(U*). 
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FIG. 2: a) Scattering amplitude for state-changing collisions, 
b) Distribution of relative kinetic energies F(E le i) for different 
values of ha, Hb, and T (in units of Eo). 

State-swapping leads to molecules in ipo,o (fifi) m 
layer B (A), see Fig. [T] For short times, the number 
density rid of these "defects" satisfies the rate equation 



drid/dt = KriAriB, 



(8) 



where ua and ns are the densities in layers A and 
B. The rate K is determined by the cross-section for 
state-changing collisions er sc , averaged over the vel ocity 
distributions Ja.b of particles in layer A and B |13j : 
K = (va sc ) = JJ v<T(p)fA{v A )fB{vB)d 2 VAd 2 VB, where 
v = v^-vg is the relative velocity for a pair of molecules 
and p = (m/2)v. The relation to the scattering ampli- 
tude in two dimensions is va = 4— f I 2 . 

m/2 l J 1 

For short times or for a small defect density, Pauli- 
blocking effects can be neglected. For reactive molecules, 
in fact, the defect density is always small because 
swapped molecules are quickly lost, as discussed below. 
For equilibrium distributions which depend only on en- 
ergy, the collision rate K becomes: 



K = 8* 

m 



|/ SC (Arel)| 2 F(E Tel )dE reh 



(9) 



where E Ie i = h 2 k 2 cl /m = ^(m/2)v 2 . Here 
g?7 



F(E lel ) = J ^-J fA(E A )f B (E B )dE cm (10) 

is the distribution of relative kinetic energies, where 
cos "f = W"' V = + v b)A e a.b = rnv AB /2, and 
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FIG. 3: Swapping rate K for different T as a function of: (a) 
the density n and (b) the imbalance /ia/mb for fiB = 0.8-Eo- 

E cm — mV 2 is the center-of-mass kinetic energy. The 
Ja,b are normalized to J °° f(E)dE = 1. 

The distribution of relative energies F(E) can be cal- 
culated from Eq. (j 1 0(1 . assuming a Fermi distribution in 
each layer f A M E ) = J\f^ 1 (e^ E ~ fJ - ( - T '>^ kBT + 1)"\ with 
the normalization constant Af = Tlog(l + e^l kBT \ The 
chemical potential n(T) is that of the ideal 2D Fermi 
gas, i.e. the solution of log(l + e ^/ kBT ) — ep/ksT for a 
given temperature and a given Fermi energy ep — fj,(0) = 
2irnh 2 /to, where n is the molecular density in the layer. 
Figure Wp shows F(E) for different temperatures. For 
temperatures fc^T 3> £f, the relative energies follow a 
Boltzmann distribution, F(E) = l/k B T exp(-E/k B T). 

The swap rate K, as given by Eq. (|9|), is the overlap in- 
tegral between F(E re \) and |/ sc (£7 re i) | 2 . The dependence 
of the collision rate on the density, inter-layer spacing 
A, and temperature can be obtained by analyzing this 
overlap. The scattering amplitude |/ sc | 2 has a maximum 
around E le i/E = (Xk re \) 2 = 0.5, see Fig. [2^, so at T = 
the maximal swapping rate occurs for Fermi momenta 
kp\ ~ 1. Hence, Fermi degenerate molecular samples 
show a peculiar non-monotonical behavior of the swap- 
ping rate as a function of density, or inter-layer spacing, 
as shown in Fig. [3^a). This effect is visible if ksT <C Eq, 
and it is smeared for higher temperatures, see Fig. [3ja). 

The swap rate K presents as well an interesting depen- 
dence on the population imbalance between both layers. 
At T — the function F(E), which may be analytically 
evaluated, depends on the chemical potentials \xa and ^b 
of the layers (or, equivalently, particle densities), acquires 
the form of a smoothed step function, see Fig. (Hb) . In- 
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FIG. 4: Swapping rate K as a, function of temperature, for 
ef/Eo = 1. The dashed line is the result of Eq. fTTj) . 

terestingly, as shown in Fig. G^b), the swap rate K may 
show a maximum as a function of the imbalance \iaI V-Bj 
if Mb i$ Eq. Again, as can be seen in Fig. GUJb), at finite 
T, the maximum persists although smeared out. 

The previous results suggest that the state swapping 
collisions may lead to very relevant effects even in the 
high-temperature regime, especially in the context of re- 
active molecules. Typical experimental values for E^/ks 
and tp I Ub are in the range of lOnK, while currently avail- 
able experimental temperatures are much higher, namely 
in the range of several lOOnK. In the following we dis- 
cuss the particularly relevant scenario of KRb molecules, 
assuming a molecular density of 5.6 x 10 7 cm~ 2 , which 
gives sf/Eq — 1.0. Figure S] shows the temperature de- 
pendence of the swapping rate under these conditions. 
For ksT ^> ep, it follows from Eq. (0 that the swapping 
rate acquires the form 



K(T) 



_ u2 \og(k B T/E ) 



k B T/E 



(11) 



Once chemically reactive molecules collide within one 
layer, they undergo a reaction of the form KRb — ► 
K2 + Rb2 and leave the trap [3] • For weak electric fields 
(no significant polarization of the dipoles), this results 
in a loss rate of K p ~ 2 X lCT 7 cm 2 /s for T = 800nK 
due to p-wave scattering^ scaling linearly with T due to 
Wigner threshold law [4j. The swapping process dis- 
cussed here leads to an additional loss mechanism for 
reactive molecules, since the creation of defect molecules 
(state ipo t o in layer B and state ipi^ in layer A, see 
Fig.[TJ allows for intra-layer s-wave collisions. Their typ- 
ical rates are K s ~ 10~ 5 cm 2 /s [3], i.e. at least one order 
of magnitude larger than the swapping rate, see Fig. |3] 
Since K s 3> K, the defect molecules are quickly lost (the 
instantaneous number of defect molecules at short times 
remains at any time very small rid — (K / K s )ha)- Hence 
we may equate the swapping rate with molecular loss 
rate, driA/dt ~ — (K + Kp)n\, assuming balanced mix- 
tures with riA — Tib ■ Note that the swapping mechanism 
discussed here becomes the dominant loss mechanism at 
low-enough T characterized by K p <C K. Since swap- 
induced losses decrease with T while losses associated to 



p-wave scattering increase linearly with T, there is a tem- 
perature threshold T c at which both loss mechanisms are 
equally relevant. From Fig. [3] and Ref. [3], we estimate 
this threshold at approximately T c ~ 530nK~ AOtp/ks 
[l5| . Therefore even in the absence of external polariza- 
tion, non-local dipole- induced swapping dominate over p- 
wave losses even deeply inside the non-degenerate regime. 

In summary, polar molecules in bilayer systems in a 
different rotational state in each layer provide a distinc- 
tive setting for the study of non-local state-swaps. We 
have shown that the special features of the inter-layer 
dipole-dipole interaction lead to a non-trivial dependence 
of the swap rate with the parameters of the system, espe- 
cially within the degenerate regime. Moreover, for chem- 
ically reactive molecules, the swapping leads to losses, 
since a swap is followed by a quick recombination due 
to an intra-layer s-wave collision. As a result of that, 
non-local state-swapping should lead to significantly en- 
hanced losses below a temperature threshold, which lies 
well inside the non-degenerate regime. Hence, our re- 
sults show that strong non-local dipolar effects should 
be observable in current experiments in the absence of a 
significant external polarization. 
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A sufficiently large electric-field gradient d£/dz may in- 
hibit state swapping. This may be neglected if Ed ^ 
2fjp giving d£/dz < 0.5kV/cm 2 for E = 0.4kV/cm. 
A Zeno-like effect, similar to that discussed by N. Syassen 
et al., Science 320, 1329 (2008), may reduce the swapping 
rate. We expect this correction to be small due to the 
moderate K s /K ratio. 



